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, We study a fundamental issue in cosmology: Whether we can rely 

\^ [ on a cosmological model to understand the real history of the Universe. 

■ This fundamental, still unresolved issue is often called the "model-fitting 
problem (or averaging problem) in cosmology" . Here we analyze this issue 
with the help of the spectral scheme prepared in the preceding studies. 

Choosing two specific spatial geometries Q and Q' that are very close 
O . to each other, we investigate explicitly the time evolution of the spectral 

distance between them; as two spatial geometries G and Q' , we choose a 
^jrj' flat 3-torus and a perturbed geometry around it, mimicking the relation of 

a "model universe" and the "real Universe" . Then we estimate the spectral 
distance between them, dN{G,Q'), and investigate its time evolution ex- 
^ ' plicitly. This analysis is done efficiently by making use of the basic results 

■ of the standard linear structure-formation theory. 
We observe that, as far as the linear approximation of the geometrical 

perturbation is valid, d]sf{Q,Q') does not increase with time prominently, 
rather it shows the tendency to decrease. This result is compatible with 
the general belief in the reliability of describing the Universe by means of 
a model, and calls for more detailed studies along the same line including 
the investigation of wider class of spacetimes and the analysis beyond the 
linear regime. 



1 Introduction 

We study in this paper one of the fundamental issues in cosmology: How pre- 
cisely a cosmological model traces the time evolution of the real Universe. This 
significant problem has been raised and studied for a long time but we have not 

^ Present address. E-mail: mseriu@eduOO.f-edu.fukui-u.ac.jp 
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yet caught a glimpse of a satisfactory understanding of this issue. It is often 
called the "model-fitting problem in cosmology" or the "averaging problem in 
cosmology" 0, 0. 

This problem originates from the very nature of cosmology. Cosmology is our 
attempt to grasp our Universe by reference to a model. The real Universe is full 
of complexity while a model, which serves as a format of perceiving reality, is 
much simpler than the reality itself. In cosmology, thus, there is in principle an 
inevitable discrepancy between the reality and its model. 

When we try to guess the past and future of the Universe based on its model, 
thus, we need a guarantee that the inevitable discrepancy between the reality and 
its model does not develop so much within the time scale of concern. Considering 
the highly nonlinear nature of the Einstein equation, the enlargement of the dis- 
crepancy would not be a very surprising result even if it were the case. Therefore 
this issue should be thoroughly clarified before we state anything meaningful for 
the past or future of our Universe. One might symbolically state the problem 
depicted here as "/s cosmology possible?' 

The central difficulty of dealing with the issue of inevitable discrepancy may 
arise from the following two facts: 

(Dl) It is difficult to define a reasonable procedure of "averaging geometry" in 
a self-consistent manner. For instance, it is difficult to define the averaged 
spatial metric (hab) in a spatial-diffeo-invariant manner. 

(D2) On one hand, we need to handle two or more universes (viz. the real one 
and its models) at the same time; on the other hand, there has been no 
efficient mathematical language for comparing two or more geometries with 
each other. 

At first sight, the difficulties (Dl) and (D2) look mutually independent. When 
we contemplate on them deeper, however, we realize that they are both rooted in 
the issue of how to quantify an intuitive statement "the two shapes look similar" : 
The issue (Dl) is reduced to the issue of how to relate the averaged geometry 
to the original geometry in a reasonable manner. The issue (D2) is about the 
comparison procedure for two given geometries. Extracting the essence from (Dl ) 
and (D2), then, we infer that the concept of "closeness" between two geometries 
would play a key role. We pursue this viewpoint further in this paper, making 
use of our preceding related studies. 

In preceding work we have indeed prepared a framework in which the mea- 
sure of closeness between geometries plays a central role P, |^, |^ , and have made 
some preliminary investigations on the model-fitting problem based on this frame- 
work 0, 0. (For simplicity we call the framework spectral scheme hereafter). The 
purpose of this paper is to make a more explicit investigation by the spectral 
scheme on the time evolution of the discrepancy between two nearby geometries, 
Q and Q', by choosing concrete, tractable geometries for Q and Q'. 
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Before going into details, let us first discuss why handling the discrepancy 
between two universes is so difficult and how we try to tackle this problem by 
introducing the spectral scheme. 

For this purpose the following example may clearly illustrate the situation. 
Let Qo be some topologically complicated space, with various small topological 
handles attached p, |^. When we observe Qq at the energy scale E, the handles 
smaller than E^^ would not be observed, resulting in a simpler effective geometry, 
say Qe- If the energy scale is decreased further to E' {E > E'), the handles smaller 
than E'~^ effectively disappear, resulting in a much simpler effective geometry 
Qe'- a similar phenomenon occurs when we change observational apparatus, 
instead of changing E; a portion of geometrical information incompatible with the 
observational apparatus would in effect disappear so that the resulting effective 
geometry depends on the apparatus used. 

In this way, we are led to the concept of effective geometry, viz. geometrical 
structure varying as a "function" of observational energy E and observational 
apparatus A, which may be indicated by G{e,a)- We intuitively know that G{e,a) 
is effectively "close" to Qq and that Q(e,a) can serve as a model for Qq when we 
are interested in the phenomena labeled by {E,A). 

On the other hand, Q{e,a) (for some E and A) is totally different from the 
original geometry Qo from the viewpoint of traditional mathematical theories of 
geometry and topology. Even though Q{e,a) and Qo are physically "close" to each 
other, mathematics classifies them into different topology classes. Topology is a 
scale-free concept and the scale of a handle is not counted, while in spacetime 
physics it should also be taken into account because the observational energy 
scale enters into the argument |]^, ^ . 

The central difficulty of analyzing the issue of discrepancy is, then, that we 
are lacking in a scheme of geometrical approximation which makes it possible to 
regard Q(e,a) as an approximation (a "model") of Qq, even though they are totally 
different in a conservative sense. Thus the first thing to do is to construct a new 
framework which justifies our intuition that Q(e,a) is "close" to in a certain 
sense. In this way we return to the issue of "closeness" discussed after (Dl) and 
(D2). 

With the above consideration, it is now clear that the following points should 
be clarified; 

(a) How to quantify "discrepancy" or "closeness" between two geometries. 

(b) How to take into account the "coarseness" of truncating information of 
geometry. (Scale-dependence of effective geometry.) 

(c) How to take into account the type of apparatus used for collecting informa- 

tion of geometry. (Apparatus-dependence of effective geometry.) 

It has been shown that a measure of closeness between two geometries, say Q 
and Q' , can be defined in terms of the spectra 0]; here the term "spectra" means 
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a sequence of eigenvalues of a certain elliptic Hermitian operator, numbered in 
an increasing order. (We mainly consider the Laplacian A as a typical Hermi- 
tian elliptic operator though the basic idea itself is universal.) The measure of 
closeness is defined by comparing the spectra for Q with those for Q' {Issue (a)). 
Furthermore, only the first N spectra (Ai, A2, ■ ■ ■ , Aat) are compared, and in this 
manner, a natural cut-off scale A^^^^ is purposefully introduced. This procedure 
can be symbolically described as "comparing the sound of Q with the sound of 
G' ". 

First of all, the spectra are the spatial-diffeomorphism invariant so that they 
are desirable quantities for quantifying "closeness" (see (Dl)). Next, on dimen- 
sional grounds, the lower (higher) spectra reflect the larger (smaller) scale prop- 
erties of the geometry. By comparing only the low-lying spectra up to the A^-th 
spectrum A^v, thus, Q and Q' are compared in a coarse-grained manner, neglect- 
ing the difference smaller than the cut-off scale corresponding to Aat. Thus the 
measure of closeness dN{G, Q') introduced in this way naturally takes care of the 
scale-dependent nature of spatial geometry {Issue (h)). Finally, the choice of 
the Hermitian operator determines the type of vibration modes used for probing 
geometry, so that in principle it corresponds to the observational apparatus used 
for measuring geometry. Thus the choice of the operator defines which aspects of 
geometry are compared to measure the closeness between the spaces. As a result, 
our measure of closeness becomes apparatus-dependent, too {Issue (c)). 

In this manner, we are now capable of quantifying the "closeness" between two 
given geometries as a "function" of the coarse-graining scale and the observational 
apparatus [Q. For a given geometry Q\ the optimal model of Q' among a set of 
models can be defined as the model geometry Q closest to Q' when measured by 
d^- Therefore we realize that a suitable model for a given "reality" Q' is not 
unique, but it varies in accordance with the observational scale and apparatus 
1^ . The mapping Q' ^ Q here may be interpreted as "averaging the geometry ^'" 
if one wishes. The unique feature of this approach is, however, it does not resort 
to the ambiguous averaging procedures such as defining the "averaged metric" 
{hab)] rather the whole of the geometry is directly mapped to some suitable model 
geometry. 

Furthermore, it has been shown that the space of all spaces Sn, which is a set 
of all compact Riemannian geometries equipped with the measure of closeness d^-, 
forms a metrizable space This result makes it justified to call the measure of 
closeness dN{G, Q') as the spectral distance between Q and Q' . At the same time, 
it means that we have constructed a basic arena Sn for analyzing fundamental 
problems in spacetime physics, such as the one discussed in the present paper. 

One of the advantages of the spectral scheme is that the basic quantities {A„} 
carry clear, well-defined meaning both physically and mathematically. Physically 
they represent the vibration modes of the space when "tapped" by a certain appa- 
ratus (A); mathematically they are the eigenvalues of a certain elliptic operator 
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defined on a space and their investigation forms one active subject in geometry. 

In reality, it is true tliat only a restricted number of cases are known for which 
whole of the spectra, {A„}^i, can be obtained explicitly. Even so, it is still mean- 
ingful to construct a general framework based on {A„}: First, let us recall that we 
do not need the whole of the spectra but only the first N spectra (Ai, A2, ■ ■ ■ , Xn) 
to be compared. It is known that the lower spectra are relatively easier to obtain 
by, for instance, the numerical calculations as well as analytical methods [ |T0(] , 
and these lower spectra contain enough information on the global nature of the 
space. Second, in many cases, we are more interested in the difference of the 
spectra {6Xn) between two geometries (or between two neighboring time-slices 
in one spacetime) than the spectra themselves, and usually we can estimate the 
former even when the latter are not explicitly available. Third, it is meaningful 
that the spectra are theoretically well-defined and in principle determined, even 
though not determined explicitly in many cases. f\ Finally, the situation is not 
proper to the spectral scheme; whatever we choose as a framework of analysis, 
investigations of complicated geometries do not go straightforwardly anyway. 

It is now possible to picture the situation in question as follows: Each point 
in the space indicates a certain geometry viewed with the apparatus A at the 
scale Aat.Q A set of models is represented by the points distributed throughout 
Sn- The optimal model Q is then represented by the point among them closest 
to the point representing the reahty Q' [|]. Thus, the real geometry Q' and its 
model geometry Q are represented by two points that are very close to each 
other in S^- These two points move in time and follow two trajectories in S^, 
and the relative behavior of these two trajectories is crucial for a model Q to 
be a good approximation of Q'. What is essential is, then, to investigate the 
time evolution of the spectral distance d^lQ, Q'). For this purpose, the evolution 
equations for the spectra of the universe is required since the spectral distance 
is defined solely in terms of the spectra. The spectral evolution equations have 
been indeed derived^, and they can be regarded as the spectral representation 
of the Einstein equation. Now we have constructed a formalism which consists of 
the following triad; (I) the spectral distance d]\f, (II) the space of all spaces 
and (III) the spectral evolution equations. 

In this paper, we apply the spectral scheme explicitly to a concrete, tractable 
situation to get deeper understanding of the relation between the reality and its 

^Indeed there are already several numerical calculations |^ (based on the finite element 
method) which provide the first 15 spectra or so for the case of 3-dimensional hyperbolic 
geometries with non-trivial topologies. This kind of numerical techniques are certain to develop 
in the near future. 

^ The situation here is very similar to the one for quantum field theory on a curved spacetime. 
There are only restricted number of spacetimes available for which we can explicitly obtain 
the complete set of positive frequency solutions of the field equation. However the general 
framework of quantum field theory is still worth constructing. 

''More rigorously, it indicates a class of geometries that look similar to each other when 
viewed with the apparatus A at the scale Xn- 



5 



model. We construct a tractable pair of spaces imitating the relation between 
the real Universe and its model, and estimate the time evolution of their spectral 
distance explicitly. As the first investigation in this direction, we here investigate 
only the linear-regime, viz. the period when the discrepancy can be regarded 
as small. Though the situation analyzed here is inevitably limited, the present 
analysis would give us insight for the cases when the discrepancy is larger than 
the present case and would provide a motivation for further studies for these 
cases. 

In Section H, we will summarize the basic formulas in the spectral scheme only 
to a necessary extent for the analysis in this paper. In Section |], we construct 
two spacetimes, mimicking the relation between the Universe and its model; as 
for their spatial sections, a flat 3-torus is chosen for what is regard as a "model 
universe" (Q), while a perturbed geometry around Q is chosen for what plays the 
role of the "real Universe" {Q')- Then we prepare several fundamental quantities 
and formulas for these spacetimes in terms of the spectral scheme. We combine 
the standard linear-perturbation theory with the spectral scheme in Section ^, and 
derive the time-dependence of various coefficients necessary for our analysis. With 
these preparations, we estimate the time-development of the spectral distance 
between Q and Q' in Section ^ (the main results are Eqs.p6| ) -(|49| ) ). Section ^ is 
devoted for several discussions. 



2 Fundamental formulas in the spectral scheme 

Let us recall some results in the spectral scheme necessary for later analysis.^ 

Let us consider a geometry Q described by a (D — l)-dimensional compact 
Riemannian manifold without boundaries, (S, h)^ 

We consider an eigenvalue problem of the Laplacian A on (S, h), Af = —A/. 
Then we get the spectra, viz. the sequence of eigenvalues arranged in an increasing 
order, {A„}„=o,i,2,--- •= {0 = < Ai < A2 < ■ ■ • < A„< ■ ■ ■}. At the same 
time, we get the set of real-valued eigenf unctions corresponding to the spectra, 
{/n}n=o,i,2,-, normalized as (fm, fn) ■= 4 fm fn y/ = Smn , where the natural 
integral measure on (E, h) is indicated by ^ := ^det{hab)- 

Now we consider two geometries Q and Q'. Let {Xm}m=o {^'n}'n=o 
spectra for Q and for Q', respectively. Then the spectral measure of closeness 
between Q and Q' of order N, d]si{Qi G'), is defined as @, || 

d^(^,^'):=E-^(f j > (1) 

^ The reader is advised to refer to Refs. Q and ||^ for more details. 
^ In this section, D indicates the dimension of a spacetime, so that its spatial section becomes 
{D — l)-dimensional. 
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where (x > 0) is any continuous function satisfying 

^(1) = , 

Tiy) > J^{x) for y > X > 1 , 
J^{l/x) = J^{x) . (2) 

For practical applications, it is convenient to assume further (i) jF(x) is smooth 
at X = 1 (then = and > from the conditions @) and (ii) 

> 0. The postulation (ii) is for making dN{G,G') sensitive enough to 
detect a fine difference between Q and Q' when they are very close to each other 
(see Eq. (|TU|)). 

It is convenient to set J-' to be J^i(x) = ^\n^{y/x + 1/y/x). Then Eq.(|l|) 
becomes § 

MO.O,^lhn0^-^) . (3) 

Now, let Riem be the space of all {D — l)-dimensional, compact Riemannian 
geometries without boundaries. We then consider a space {Riem, d^)/^, viz. the 
space Riem equipped with the measure of closeness dN (Eq-dl)) along with a 
natural identification (~) of isospectral manifolds.[| 

It is proved that the space of all geometries {Riem, dN)/^ forms a metrizable 
space @].| It means that the measure of closeness dN (Eq.(||)) can be regarded as 
if it were a distance in spite of its mild violation of the triangle inequality 0, ^. 
From now on, thus, let us call diq{Q, Q') in Eq.(^ the spectral distance of order 
between two geometries Q and Q' . Let Sjq be the completion of {Riem, djy)/^ by 
means of the latter's metrizable structure, since such a complete space is math- 
ematically more desirable. One might call Sn the space of all spaces of order 
N, or the spectral space of geometries of order A^. Mainly due to its metriz- 
able nature, Sn possesses several nice properties such as the second countability, 
paracompactness (then the partition of unity can be introduced on S^) and 
locally-compactness (then an integral over Sn can be defined) .0 

When we focus on the discrepancy of two geometries that are very close to 
each other in Sn, there is no significant difference even if we replace J-'i{x) (used 
in Eq.(^) with some other continuous function J-'{x) satisfying the conditions in 
(ID (see Eq.(|l^) below and the arguments there on this point). 

Let us introduce some convenient notations for later use. Let A and Aab be 
any function and any symmetric tensor field, respectively, on a spatial geometry 

^When two Riemannian manifolds possess the identical spectra of the Laplacian even though 
they are not isometric, they are said to be isospectral to each other. There are some known 
examples of the isospectral manifolds. For more details, see Refs. [|lO| and For the physical 
interpretation of isospectral manifolds, see Refs. § and §. 

^ It turns out that [Q the distance function for metrization is provided by Eq.(^ with the 
choice !Fo{x) :— ^ \nmax{y/x, l/y/x) for !F. 

^ For basic facts of point set topology, see e.g. Refs.[^ and p3[ . 
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(E,/i). Then, we define spatial-diffeomorphism invariant quantities (A)^ and 
Ua&>„ as 

{A)„ := / fnAfn^, {Aa,)„ := ^ f d'^ fn Aa, U V ' 

where is the n-th eigenfunction. For a quantity (w4a6)„, it is assumed n > 1. 

Any function A{-) on the spatial section E can be expanded in terms of 
{fn}n=o^ such as A{-) = E^=o A /n(-) ■ Here we note 

An = ifn , A)-= [ fn A^ . (4) 

The 0-th component is related to the spatial average of A over E, 

where V is the {D — 1)- volume of the spatial section (E, /i). 
We also introduce the quantity 

(/ m n) := = y" , (5) 

which is totally symmetric in Z, m and n. 

It is useful to define eat and rab, the trace-free components of the extrinsic 
curvature Kab and the Ricci tensor Rab for (E, h) , respectively; 

1 1 

^ab ■= - Khab , Tab ■= Rab " ^ _ ^ Rhab ■ (6) 

These quantities describe anisotropy of the geometry (E, h). 

At this stage it is appropriate to make it clear what the terms "the real 
Universe" ("reality") and "the model" actually indicate in this paper. We shall 
compare two spaces, Q and Q', that are very close to each other in Sn', this is 
a totally well-defined situation without any ambiguity. Throughout this paper, 
we deal with only this mathematically well-posed problem. Only on the final 
stage, however, we may interpret the results of this analysis in the context of the 
averaging problem in cosmology, supposing that one of the spaces corresponds 
to the real Universe and the other to its model. To keep in mind this final 
interpretation, we symbolically refer to the simpler space among the Q and Q' as 
the "model" , while the other space as the "real Universe" ( "reahty" ) . 

It is now in principle possible to investigate the evolution of the discrepancy 
between the reality and its model generally. In this paper, however, we restrict 
ourselves to the cases when 

(Al) the discrepancy between Q and Q' can be totally described in terms of the 
difference in the spatial metric, and 
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(A2) they are initially so close to each other in Sn that the linear treatment of 
the discrepancy is justified at least during a certain period of time. 

(The assumption (Al) implies that the two spaces possess the identical global 
topology.) 

Now, let us consider two nearby spaces Q and Q' in Sn] their discrepancy is 
represented by a small difference in their spatial metrics, 

7afe := Kb - hab , (7) 

where hab and /i'^^ are the spatial metrics for Q and Q', respectively. We treat •jab 
as a small quantity and indicate its order of magnitude by 0(7). 
By means of the variation formulas for the spectra |p, we get 

^lnA.:=^^ = -(7a.)„-^(7)„ , (8) 

where 7 := /i^Vafe and j^f, := Jab - \lhab. 

Since Eq.(|D plays the role of our key equation, let us briefly review how to 
derive it. The discrepancy of metrics, -^ab '■= Kb " ^ab, causes the variation of 
the Laplacian, 6A, which induces the discrepancy of the n-th spectrum, 5A„ : = 
A'„ — A„,. First we recall there is a well-known relation ("Fermi's golden rule") 

<5A„ = -(5A)„ := - fJMn^ . (9) 

Thus we need the expression for (5A. For this purpose we take the variation of 
both sides of A/ = ^~^9„( h"-^dhf\j) with respect to hab, yielding 

5A/ = Uald'^f - l-daiYcd-fV) ■ 

^ V 

Thus we get, 

(5A)„ = 1| f^daicffn^- J fnda{l\d^fn^) • 

By suitable partial integrals, this expression along with Eq.(D yields Eq.(||). 

Let us now consider the general form, Eq.(|l|). We insert A'„ = A„ + 5\n into 
the R.H.S. (right-hand side) and expand it in terms of 5A„. With the help of 
Eq.(H), then, the leading term in Eq.(|]) becomes p 

ciiv(^,^') = ^-^"(l)E((7j. + ^(7)„)' • (10) 

We thus observe a prominent feature; whenever Q and Q' are very dose to each 
other in Sn, the leading term of the spectral distance dN{Q,G') is universally 
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given by Eq. (]7^, irrespective either of the detailed form of the spectral distance 
or of the gravity theory. 

We note that the choice of the function J^{x) affects only the unimportant 
numerical coefficient of d^iG^Q') in the leading order, irrespective of the cut-off 
order A^. Thus, it suffices to consider a particular form given by Eq.(|^). Then 
= |, so that Eq.dig) becomes 

1 ^ / 1 \2 1 

where in the last line, 7 indicates a vector in whose n-th component is 
In '■= {lab)n + |(7)n- further derive the expressions for d]^]{Q,Q') and 

dN{Q,Q') (For more details, see Ref.0). 

Equation Eq.(p!lD plays the role of the central equation for our investigation 
on the evolution of the discrepancy between Q and Q'. 



3 Quantitative description of the relation be- 
tween two nearby geometries 

3.1 Models to be investigated 

Let us now investigate explicitly the discrepancy of two specific nearby geome- 
tries.H 

As a "model" universe, we construct a spacetime in the form of x R as 
follows]^: Let (a = 1,2,3) be the Cartesian spatial coordinates on R^. By 
imposing the identification ~ + 1 (a = 1, 2, 3), then, x" (a = 1, 2, 3) turns 
to a coordinate on a circle 5*^, taking its value in [0, 1] with the identification 
~ 1. We then consider a metric defined by 

ds"^ = -dt^ + a{tf5abdx''dx^ , (12) 

where a{t) is the scale-factor of the model. The t = constant spatial section of 
this model, S^, corresponds to Q in the previous section. 

As for the "real Universe" , a perturbed geometry around the above model is 
chosen: We prepare a spacetime x R with the same spatial coordinates as 
above, x'" G [0, 1] with ~ 1. The metric is then given by 

ds"^ = —dt^ + habdx'°'dx'^ , (13) 

We set D — 4 hereafter unless otherwise stated. 
^^In this paper, Roman indices (such as a, 6, • • •, k, I, ■ ■ •) indicate the spatial indices and 
they are raised and lowered by the spatial metric. On the other hand, Greek indices (such as a, 
/3, ■ ■ •) indicate the spacetime indices and they are raised and lowered by the spacetime metric. 
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where the spatial metric hat is given by 

hab{t,x"') = a^{t) hab{t,x"') := a^{t) {6ab + ^ab{t,x"')'j , 

7a6(t,a;'"): = a\t)%bit,x"') . (14) 

The t = constant spatial section of this spacetime, S'^, corresponds to Q' in the 
previous section. 

For later convenience, let T^(l) denote the regular 3-torus that is confor- 
mally equivalent to Q; T^(l) := [0, 1] x [0, 1] x [0, l]/~, with the standard metric 
dP = 5abdx°'dx^ (/^ indicates the identification for making the regular 3-torus). 
Let {A„}„=o,i,2,-- and {/n}n=o,i,2,--- are, respectively, the spectra and the eigen- 
functions of the Laplacian on T^(l), A := d'^^dadh- 



3.2 Spectral representation of geometrical quantities 

In view of Eqs.([A7|), ( |A11| ) and ( |A16| ), the metric perturbation 7^6 in Eq.(|l^) can 
be expanded in terms of fAhab, 5"^, Cat w^^: 

%b = \T. h Kb + E ^^b + T:r2 ~Cb + E ^ab ■ (15) 

A A A A 

Here the prime symbol in J2a implies omitting the zero-mode of A from the 
summation; the coefficients 7'^'', 7'^', 7'^^ and 7'^' are attributed to the pertur- 
bations of conformal-type, scalar-type, vector-type and tensor-type, respectively. 
We consider any perturbed quantity up to 0(7) hereafter. From Eq.(|T5|), we get 

7 = /^"Sa. = E7'A' h . (16) 



Now from Eqs.(|T4D and (|15D , it is straightforward to compute Rab, R and Tab 
(see Eq.(^) for Q'. Extracting the /^-components of R and D^-D^Tab (according 
to Eq.(l)), thus, we get 

Ra = IxAiY'^ + (A^O), 

{D'^D'T^b)A = -^Aa(7<^'+7'^')v^ (^7^0). (17) 

The two formulas in (|17[) are equivalent to each other due to the Bianchi identity, 
D^iRab - \Rhab) = O.g 

In the same manner, K and eab (see Eq.(|])) are expanded as 

K = K^,{l + Y^KAfA) , 
A 

eab = ^av(E'(«'/5r)^f. + E'(«W)CUE(«'A 
{ A A A 



^"^ For, the spectral representation of the Bianchi identity becomes Ra + ■^{D"^D^Tab)A = 
(A^0)@. 
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Extracting the /^-components of K and D°-D^eab, thus, we get 

^ r 3HVV iA = 0) , 

iD''D''eab)A = 2H\A(3fy/^ (A ^ 0) . (19) 

Here we used the relation i^av = 3if := 3a/a; To obtain the last equation, the 
second equation in (|A8|) has been used. 



3.3 Spectral representation of matter quantities 

We consider the perfect-fluid type matter, T"^ = (p + p)u°'u^ + pg'^^, where 
u"" is the 4- velocity of the fluid element {u°'Ua = —1). The matter velocity, u°' 
(a = 1,2,3), is regarded as 0(7). The standard parameters are introduced for 



later use: u := p/ p {0 < v < 1/3 for normal matter) and Cs := \J dp/dp, the 
sound velocity of the matter relative to the fluid element. 
We expand p and p as 

P = Pav(l + 5), P = Pav + PavC^f^ , S = Y,'SAfA ■ (20) 

A 

Then, the spectral representation of p becomes 



~ I 5^PavV^ (A ^ 0) . ^''^ 

The expansion of matter 6 := DaU°- is also expanded as 6' = SHJ^A^AfA, then 
its spectral representation becomes 

eA = SHeAV^ (A^O) . (22) 

Then the divergence of the momentum density of matter, = {pav + Pav)ua + 
0(7^), reduces to the spectral representation, 

{D ■ J) A ^ (Pav + P..)0 = 3(1 + iy)Hp,,QA^ (A^O) . (23) 

3.4 Spectral representation of constraint equations 

The Hamiltonian constraint is represented in the spectral form as 

Ra + IT. Ka'Ka"{A' A" A) - -PA - 2KVV5ao - (eafce''')^ = , (24) 

O A' A" ^ 

where a := (See Eq.(|D for the definition of {A' A" A).) We neglect the 

quantities of order 0(7^). Then, for A = 0, Eq.(^) becomes 

6ij2 - -Pav - 2A = , (25) 
a 
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with the help of Eqs.(|T|) and (Here we note {A' A" 0) = -^6a'A"-) Need- 

less to say, Eq.(p5D coincides with the standard equation for the (locally) flat 
Friedmann- Robertson- Walker model. On the other hand, Eq.(p^ becomes, for 

^A^(7r+7r) + 12//V--Pav5A = (A^O) , (26) 
6 a 

with the help of Eqs.([l3), (0) and (|2l|). (Here we note Ea'A" Ka'Ka"{A' A" A) = 

ISH'^kaVo? + 0{j^) for A ^ 0.) With the help of Eq.(^), Eq.(^) is represented 

as 

;)6a-2ka\ (A^O) . 



(C) I (S) 

Ta +rA 



A. 



(27) 



On the other hand, the momentum constraint in the spectral form becomes 



Ka 



3 1 
2X2 



{D'^D''ea,)A + —{D-J)A\ = Q (A^O) 



With the help of Eqs.(|l9]) and (|2|), thus, Eq.(|D gives 

{ka + PT) + + ^)(3i^' - A)e^ = (A ^ 0) 

Z Aa 



(28) 



(29) 



where Eq.(p5D was used to eliminate pa 

The coefficients k^'s and /3^'s are related to 7^'s: From Eqs. (|13D and (|l^ . 



we get Kab = \hab =f hab + \l ab- Together with Eqs. (|T5|) and (|T6[), this relation 
yields the spectral expressions for K and up to 0(7). Comparing them with 
Eq. (|TB|) , we finally get 



a 



--: 66H 



(S) 



6H 



(A^O) 



6H 



iT (^7^0) 



6H 



Then Eq. (|29|) reduces to 



9H 



ilT+lT) + —{l + iy){3H'-A)eA = (A^O) 



7r-(30) 



(31) 



In particular, we note that 7^' + 7^' = constant when = 0. 

We have obtained necessary relations to investigate 7^ := (7ab)^ + a^ 
the A-th component of 7, which is important in view of Eq.(|l^). With the help 
of Eqs.(|l|), ([T|), (^, (^), (IM2D and (|7n7|) , we get after some manipulations, 

lA = 1^0^' + + + 1%') + E ^S'7y , (32) 



where c)^ji/, c^j^/ and c)^^/ are given in Eqs. (|A5|) and (|A17|) . We note that the 
vector components do not appear in Eq.(^); it can be traced back to the spatial- 
diffeomorphism invariance of the spectral distance (Eq.(|l])). Thus we can safely 
omit the vector quantities (those with suffix (V^)) hereafter. 
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4 Combination of the spectral scheme and the 
hnear structure-formation theory 

Looking at Eq.(^), the behavior of coefficients 7^', 7^' and 7^^ determine the 
evolution of the spectral distance d^iG, G')- Considering Eq. ([T5|) , these behaviors 
are translated from the time evolution of 7ab, which can be estimated by the 
structure-formation theory. (The basic results of the standard linear structure- 
formation theory is summarized in Appendix B Jl^-) Here two scales control 



the situation: The physical wave-length of the perturbation (of geometry and 
matter), Aphys, and the causal scale at the time t, Cst (c^ is the sound velocity of 
matter). 

4.1 Behavior of coefficients ^Ai Pa and i^a- The case of 
Aphys > Cgi and 6' 7^ 

We ffist consider the case (I) Aphys ^ Cgt and 6^ 7^ 0. 

Since > t in this perturbation mode satisfying this condition 

does not oscillate in effect during the cosmological time-scale t. Comparing the 
expansions of 7, 5 and 6 in Eqs.([T6|), ( pO]) and (p^) with their behavior in Eq. ( [Bl| ), 
we see that 

. .rw^-^.. , e.w^il^^.. . ,33) 

where we have used H(t) = ^.{i+u) + ^(7) (^he standard result derived from 
Eq. (|25|) ) in the last equation. 

On the other hand, comparing the expansion of 7^/ in Eq.([T5|) with their 
behavior in Eq. (p2D , we get 



1 f a \ , If Anhvs \ ^ /.N 1 



^Tit), lT{t) oc ^ (^^j 6At) oc ^ i^^j 6At) oc —t-TT. . (34) 
We omit vector quantities since they do not contribute to the behavior of the 

An 



spectral distance. Since ^ 1, Eqs.(^) and (Q) indicate 

17^1 « 17^1, 17^1 , 

which is compatible with Eq. (p3| ). 
Now from Eq.(|30D we see that 

^a = -^¥a ^i'a ^-t^ ' (35) 
where we have used [7^^ | oc |7a^^ |/t since t <^t (see Appendix ^ A\ ) and oc 1/t. 
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In the same manner, from Eq.(|30| 

^(S) ^ ^/s) ^ ^ / ^phys 
1 



V \ c.t vk^ 



oc 7j' oc -^t-TTi; . (36) 
Comparing these coefficients with 6a we see that 

i\5A\ ~ |7ri ~ « i\lT\ ~ l/^n, 17^1 ~ l/^ri) • (37) 

4.2 Behavior of coefficients 7^, and The case of 

Aphys ^ Cst and 6* = 

We next consider the case (!') Aphys ^ c^t and 6* = 0. 

Comparing the expansions of 7 and 6 in Eqs.([l6D and (|20|) with their behavior 
in Eg. (P^) , we see that 



2(1+3^) 2 



Next, comparing the expansion of 'jki in Eq.(|T5|) with their behavior in Eq. 
we get 

iTit), iTit) oc 6A{t) cx (^)^A(t) cx . (39) 

From Eqs. (pHf) and (0), we observe that 

17^1 « 17^1, 17^1 , 

which is compatible with Eq. (P7|) 

In the same manner as in § [4.1| , we get from Eq. (|30D that 

'^A = TTFjlA OC oc t 3(1+^) , 

= -^if oc u6a cc ut 3(1+^) , = — cx Pt 3(1+^) , (40) 

where we note H-^^f,H-^^^P oc v5a due to Eq.(|B6D. Comparing Eqs.(||)-(|4g) 
with each other, we thus get 

m - i7ri) - - Wa'\) - (1/3^1-^17^1), mr\-t\^T\) 

« 17^1,17^1 • (41) 
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4.3 Behavior of coefficients 7^, Sa, Pa and ka- The case of 

Aphys ^ Cst 

Finally, we consider the case (II) Aphys -C Cgt. A perturbation mode satisfying 
this condition oscillates many times within the cosmological time-scale. 

Comparing the expansions of 7, (5 and 6 in Eqs. (|T6[) , ( ]20|) and ( ]22|) with their 
behavior in Eq.( [B8|) , we observe that 

1 1-3 

6A{t) (X -^rWT^-OsC , 7^^>(t) ~ -(5A(t) , 

QA{t) oc - ^r^l^^ ■ Osc' ~ v^t^ ■ Osc' . (42) 

Here Osc and Osc' indicate the oscillatory factors with the period r = 2'Ka/kcs 
whose relative phase- difference is 7r/2; we have used H oc 1/t and assumed Cg ~ 
constant to reach the final estimation in the third relation. 

We now compare Eqs.(^) and (^) with Eq.([B9D. For instance, we can es- 
timate as 1/5^' I oc |iJ^"^7^'| oc 7|t7^'|~ (-ff^) ^ I^aI- In this manner, we 

{\Sa\ - 17^1) - ihTl » iM - ipn \PT\) ■ m 

5 Evolution of the discrepancy between two nearby 
geometries 

5.1 Some useful propositions 

Let us prepare propositions that become useful soon. 

We first introduce some notations. Let V = (Vi, V2, V3) and W = (Wi, W2, W3) 
be vectors, whose components are non-negative integers, viz. V,W e No^ (where 
No := {0} UN). 

Now, i^{V) is defined as the number of zero elements in V.f^ 

Next, we say that V and W are compatible with each other and write V ^comp 
W when jj^{V) = i^{W) and the positions of their zero elements are identical.^ 

Then we define V^\(\y) to be the component of V that is compatible with W, 

viz. Vj|(^) := {il-5w,,o)Vu il~5w,,o)V2, (1 - ^^^3,0)^3) .Q 

Similarly, we define V^_l(i4/) to be the component of V that is causing the 
incompatibility with W, viz. V|,(^) := V - V||(^) = {6wuoVi,Sw2,oV2,Sw3,oV3). 
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For example, #(1, 2, 3) = 0, #(3, 0, 2) = 1 and #(0, 0, 0) = 3. 
14 For example, (1, 3, 0) ^comp (2, 1, 0), (2, 0, 3) ^comp (3, 0, 3); (2, 3, 0) T^comp 
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For example, when V = (2, 3, 1) and W = (0, 1, 2), V||(^^ (0, 3, 1). 
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We may simply write V\\ and V± instead of Vji^^-j and Vj_i^^-j, respectively, 
whenever obvious. 

As is explained in Appendix A, the basic modes of A on the model space E are 
labeled by the pair of two vectors, {n, a), (tt-i, n2, e Nq; ui, (T2, cxs e {0, 1}). 
We often write simply A for {n,a). Note that A = (0, a) automatically implies 

A = (0, 0) since fjl^'^'^ = is the only possible eigenfunction for n = (see 
Eq.(^). 

For A = {n, a) and A' = {n', a'), the quantity # (o^||(n')) (the number of zero 
components in (^\\(n')) has just been defined above. 

Now we can state the proposition. Let caa' be any function of A = {n, a) 
and A' = (n',a') containing the factor (— Then, for which is any 
function of A, it follows that 

Proposition 1 

A [ "oEa Cao {n =0) . 

(Here and cao are abbreviations for aj-^jg) and (-(j^j-), respectively; the 

prime symbol in X)a' indicates A = (0, 0) is not included in the summation.) 
For the proof, we note that the L.H.S. (left-hand side) contains a factor 

J2a{~)'^^'^^^^"'^^ ■ This factor is nothing but the summation of -|-l's and the 
same number of — I's (so it vanishes), unless "o^||(n') = independently 
of ct" ; it means that the L.H.S. vanishes unless n' = 0, so that the 
proposition follows. 

As a corollary, it immediately follows that 

Corollary of Proposition 1 

A A' 

Now, we show another proposition. Let A = {fl, a) and A' = (n', a') once more. It 
may be clear now what riiifn') indicates, and so does 5^/ 2n„,^i^- Then this time, let 

1 1 V y 'II (.^ ) 

caa' and c'aa' be any functions which contain the factor (— )*^'^i""''''5ri',2n||(-,)5^/ o- 
Then, for oa' and /Sa' that are any functions of A', it follows 



Proposition 2 



J2'cAA'aA' Y^Caa"(^A" = Y^Y^CaA'CaA'^^A'I^A' 
A \ A' / \a" I A A' 



For illustration, suppose n' = (1,0,0). Then, <?||(n') = (fi,0,0) so that # (f?||(n')) = 
3,2 for CTi = 0,1, respeetively; then (— )*(°"ii("')) = —1,1 for cti = 0,1, respeetively. Thus 

^^,(_)#(-,K.',) = 0. 
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We note that the L.H.S. contains a factor u := Eal-)**^^"^"'^^ ("")). 
Every term in m is 1 when n' ^comp so that u becomes some (non-zero) 
natural number when n' ^comp n"] on the other hand, when n' i^comp ^"5 u is 
the summation of +l's and the same number of — I's so that it vanishes. It 
means that the factor u can be regarded as oc b^' .^n because of the presence 
of ^n',2n||,-,^'^n",2n||,-,,^ ou the L.H.S. Thus, uotiug the presence of another 

II } II ) 

factor bgi ^6^,, g on the L.H.S., the claim follows. 
Finally, suppose caa' satisfies the same condition as in Proposition 2. Then 
Proposition 3 

A A' A A A' 

It immediately follows from Proposition 1, its Corollary and Proposition 2. 



5.2 Estimating the time evolution of (1^(0, Q') 

Now, let us look at Eqs.(|llD and (^). With the help of Proposition 2 and 
Proposition 3, we can simplify the expression for 7^^. We assume that the 
coefficients 7^', and 7^' are invariant under the permutation r of A, r : 
A h-i> t{A). (More explicitly, r : {n,a) t—>- (^r(n), r(cr)^ .) Then, we note that any 
term in 7 ■ 7 of the form 

V^'V^' (T) (T) 

A A' 

vanishes unless caa' = Caa' because c\^, is odd for the permutation a : AA' h-* 
a{A)a{A') (see Eq.( |A18|) ), while other coefficients c^^\ c^^^ and c^^^ as well as 7^', 
7^\ and 7^' are even under the same permutation (see Eg. ( |A5|) ) . In a similar 
argument of permutations, the contribution of c^^q to ^ a is only in the form of 

Thus, we get 

dN{G,G') = ^7-7 , 
lo 

A A' 

Here we note that N is the number of eigenvalues less than Ajv; the contribution 

of the term J2'aCao 7o is included in the last term on the R.H.S. 
We pay attention to the three typical types of perturbations: 
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(I) The case of Aphys ^ Cgt and 9^0 
(r) The case of Aphys ^ Cgt and 9 = 

(II) The case of Aphys ^ Cgt 

For the case (I), the behavior of d^i^Q, Q') is controlled by the behavior of (7^' + 
7^,')^ ^ 7^,'^ and 7^'^, due to Eq.(^). The situation is same for the case 
(r), too, due to Eg. (p|) . On the other hand, for the case (II), the behavior 
of dN{G,Q') is determined by the behavior of all 70^^'^, (7^' + 7^'')^ and 7^'^ 
because of Eq. (^31) . 

These observations enable us to estimate the time-dependence of 7^^. First, 
looking at the explicit expression for c^a''s (Eqs.(|X5|) and ( |A18| )), the summations 
over a and a' on the R.H.S. of Eq.(|4^) yield no essential contributions to d(^Q, Q') 
within the present order of accuracy. Thus, we can safely replace A with the 
coordinate wave-length k defined as /c = 2™. Then we may write 7^ instead 
of 7a, so that dN{G,Q') is estimated as dN{G,G') ~ / •yldfi^k), where fi{k) is a 
suitable measure for k. Then the estimation of 7^^ for the cases (I), (F) and (II) 
can be summarized as (see Eqs.(|3^), (|39| ) and (^2])) 



mt-'-^ tor (I) 



7,= -^ BMt" for (I') , (45) 



fc4 

2(l-3i/) 

mrm+uT . Osc^ for (II) , 

where D{k) denotes the extra fc-dependence coming from (a) the coefficients 

i'^AA' + '^AA'Y '^AA' ; (^) ^^e spectral property of the perturbation 'jab, 
i.e. the explicit value of 7a's in Eq.(|15[). 

The fc-dependence coming from (a) is determined without ambiguity by the 
definition of caa''s (Eqs.(p^ and ( [A18| )). As is explained after Eqs.(|A5|) and 
( |A1^ ), the coefficients caa''^: do not depend on the magnitude k, but only depend 
on the relative direction between k and k' The summation of both k and 
k' {A and A') in Eq . (|4^) ensures the summation for all the relative directions k 
and k', which yields only minor corrections. The fc-dependence coming from (6) 
depends on the initial condition for ^ab- As one natural assumption for the initial 
condition for 7^^, we here postulate that the ^-dependence coming from (6) shows 
no peak at a particular scale but a mild fc-dependence. 

In any case D{k) is not essential for the t-dependence of dN{G,G') (see ar- 
guments after Eqs. (|48|) and (|4^)). We here take care of D{k) by postulating a 
power-law, D{k) oc k~^ for (I) and (F), and oc for (II). 

We can estimate diy{Q,Q') by converting the summation J2a' i^i Eq.(p^ to 
fc-integral from k = 1 (because of the structure of the space) up to /c = /cat, the 



^'^ Note that k = 27rn and k'^ = A.ir'^n ■ n (Eq.( Al) and below) 
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maximum value for k determined by the cut-off order A^. (It is a standard way 
of estimating a summation in statistical physics, and considering Eg . (|A2|) , this is 
a reasonable estimation when the universes are sufficiently large.) We consider 
two situations, when we are studying 

(A) relatively shorter scale geometry as well as the global features of the Uni- 
verse, and 

(B) only the most global features of the Universe, 

where the horizon scale at the time is the natural standard for "s/iort" or "/on(yf" . 

Situation (A): This situation corresponds to > kn- Here kn '■= 27Ta/cst oc 

t 3(1+") , the wave-number corresponding to the horizon scale at t. From 
/c = 1 to /c = /c// in the integral region, we can use the result for the case 
(I) (or (r)); from k = k^ to k = k^, we can use the result for the case (II). 
Thus we estimate, based on the case (I) for k G [1, kn], 

2(1-1-) fkn 1 „ . 2(l-3t-) fkN 1 „ „ 

dN{g,Q) ^ A't-^+^ —-k^dk + B't-^(^ -—■Osc^Pdk 

J 1 k J kjj k ^ 

~ A^{t) + BkN^-'^ 7]{t) r WIT . (46) 

Here, A and B are appropriate positive numerical factors that make the ex- 

(l+p)(l+3i^) (2-9)(l+3i/) 

pression dimension-free; C,(t) := 1 — Ct i+i^ and?7(t) := 1—C't 3(1+1^) 
with C and C being appropriate positive numerical factors that make ^{t) 
and ri{t) are positive for t > 3to- We note that ^(t) — and ?7(t) — 1 as 

2(l-3i/) 

t — > 00. Thus as time goes on, the second term with the behavior ~ t 3(1+^) 
becomes dominant. 

In the same manner, based on the case (I') for k G [1, kn], we estimate 

d^iG, G') ^ A^t) + 5A;^2-, ^^^^ ^-'-^ (47) 

Situation (B)\ This situation corresponds to /cat is less than kij. Thus, for the 
case (I), 

2(1-1^) r^N 1 „ 
dN{G,G') ^ A't-^ -^.k^dk 

^ A{1 - kj^-^^+P^) . (48) 

In the same manner, for the case (I'), we get 

dj,iG,G')^A'il-kr,-^'^^^) . (49) 
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Looking at these results, Eqs.(P6D-(^9D, we see that dN{Q,Q') decreases with 
time (with Eg . (|49|) being a marginal case of constancy). Typical behavior is 

2(l-'iv) 2(1-1/) 

dN{QiG') ~ t 3(1+1/) for the situation (A), and ~ t 1+" for the situation (B). 
It means that, for the situation (B), d^iG^Q') decreases more rapidly than the 
situation (A). This might be interpreted as that, within the linear regime, the 
dynamics of the most global structures of the Universe are approximated by those 
of a model quite well, compared to the shorter-scale features. 

Significantly, Eqs.(^6D-(^9D are also showing the scale-dependence of d^iQ, Q') 
as an increasing function of the cut-off scale k^- It indicates a reasonable result 
that the finer structures are observed, the more discrepancy between the reality 
and its model is detected. This result is consistent with our observation in the 
previous paragraph. 

Finally, we realize that D{k) has no significant influence on the t-dependence 
of di\f{Q,Q'); rather D{k) for the integral-region (II) determines the details on 
how d]^{Q, Q') depends on the cut-off scale N . 

We observed that, in the regime where the linear approximation of geometry 
is valid, d]^{Q,Q') does not increase with time, rather it shows the tendency to 
decrease. This result is favorable to the reliability of describing the Universe 
approximately by means of a model, at least during some period of time. Needless 
to say, the analysis here is limited to just one pair of geometries within the linear- 
regime. The present result calls for more investigations along the same line on 
the behavior beyond the liner-regime and wider class of spacetimes. 

6 Summary and Discussion 

In this paper, we have studied a fundamental issue rooted in the foundation of 
cosmology: To what extent we can rely on a model to recognize the Universe. To 
analyze this issue, we argued that the concept of "closeness" or "discrepancy" be- 
tween two geometries was indispensable, so that we have resorted to the spectral 
scheme developed in the previous studies. Then we have estimated the discrep- 
ancy measured by the spectral distance between two spatial geometries, Q and 
Q' , that are very close to each other, mimicking the relation of the Universe and 
its model. 

We have set ^ to be a fiat geometry constructed from the fiat Friedmann- 
Robertson- Walker model obtained by the point-identification; as for Q' , we have 
chosen the geometry generated by perturbing the geometry Q. These Q and Q' 
were regarded as the imitation of the relation of the reality and its model. 

Then the behavior of the the spectral distance d^iG, Q') has been investigated 
explicitly. With the help of the linear structure-formation theory, the behavior of 
dN{Q,G') has been estimated as Eqs.(^)-(P5|). In particular, it is an important 
observation that the geometry difference contributes to the discrepancy of the 
two universes through a particular combination 7„ := {^ab)n + (Eq-OTT|))- 
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The spectral scheme served as a powerful tool for dealing with the present 
problem. It is effective to analyze the type of problems involving two or more 
geometries. Indeed we were able to derive several non-trivial relations describing 
the discrepancy between the reality and its model which would have been difficult 
to find out without the spectral scheme. Here we realize primary importance of 
preparing a suitable language for talking about an entangled situation before 
tackling it. 

We have observed that, as far as the linear approximation of geometry is valid, 
dN{Q,Q') does not increase with time prominently, rather it shows the tendency 
to decrease. This result is, at least, compatible with our belief in the vahdity of 
the approximate description of the Universe in terms of a model. The present 
investigation is on just one pair of geometries within the linear-regime; it is now 
required to investigate wider class of spacetimes as well as to study the behavior 
beyond the liner-regime along this line. 

We have learned in this paper that it is practically possible to proceed in 
this direction; now we can improve the order of accuracy of approximations, e.g. 
the passage from Eg. (^51) to Eqs.(|46|)-(^9D. One drawback of the analysis linked 
with the linear structure-formation theory is that detailed features of the problem 
become, as it were, hidden behind it. Therefore, it is desirable to analyze the same 
problem via another, more basic route, i.e. by directly using the fundamental 
equations for dN, d^, dN etc. [Q. 

As the next step of investigation, the evolution of d^lQ, Q') should be studied 
beyond the linear regime of the geometrical discrepancy. It can be done since 
we have fundamental relations Eqs.(^) along with the spectral evolution equa- 
tions 0. 

Finally we note that our approach to the averaging problem in cosmology is 
quite different from the "usual" one which the term "averaging" might suggest. 
In the spectral scheme, there is no ambiguous procedure of "averaging geometry" 
such as (hab) (recall the difficulty (Dl) in §|I|). Rather, we first define the con- 
cept of "closeness" between geometries and choose a model Q which is a simple 
geometry and reasonably close to the reality Q' in the space of all spaces Sn- 
This mapping from Q' to Q is what we regard "averaging" in effect. Needless to 
say, this procedure is mathematically well-defined. Since the concept of "close- 
ness" (djy) is scale-dependent and apparatus-dependent, this mapping procedure 
( "averaging" ) is also depends on both the observational scale we are interested in 
and the apparatus we rely on. These are desirable properties and we can expect 
that the framework outlined here is glimpsing a new way of viewing spaces as 
scale-dependent and apparatus-dependent objects. 

The type of analysis developed here has a general applicability to a wider 
range of theories, not restricted to cosmology. It is an issue of the internal 
relation of the triad {Reality, Model, Dynamics), which arises quite universally. 
When a theory deals with objects characterized by a sequence of real numbers 
ordered in an increasing order ("spectra" of the theory), {A„}„=o,i,2,--) can 
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construct the spectral measure of closeness similar to Eq.(|I]) along with Eq.(^. 
Then one can construct the spectral space S, which is similar to our Sn- This 
space forms a metrizable space when the spectral measure of closeness (Eq.(0)) is 
chosen suitably, e.g. like Eq.(|^).0 A "good" theory should have such a property 
that two generic points close to each other in S remain close during some period 
of time of concern. Then we are entitled to approximate the reality by a model 
or classify the objects into types. The analysis presented here is the very first 
step for understanding the issue of types in physics. It is interesting to analyze 
fundamental theories from this viewpoint. 
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APPENDIX 



A The basic modes 

Here we present the basic modes on the t = constant spatial section, E^, of the 
"model" universe x R with Eq.(|T^. We recall that we adopt the Cartesian- 
type spatial coordinates (a = 1, 2, 3) which take the value in [0, 1] with the 
identification ~ 1. We note A = ^dadbS""^. 
Now, let 

_ f cos^ (f^ = 0) , 
" 1 sin^ (a = l) . 



COSn 



Then the eigenf unctions of A on Sj can be represented as 

^2(3-#('?)) 



fj[\xi,X2,X3) = I— cos^, 27mia;i cos^^ 27rn2X2 cos^3 2^3X3 

(ni,n2,n3 = 0, 1,2, ■ ■ ■; (Ti,a2,cr3 = 0, 1) , (Al) 



where a := (o"i, (72, CT3), n := (ni,n2,n3) and denotes the number of zero 

elements in The corresponding spectra can be represented as 

(27r)2^ ^ (27r)2 .... , 
XH = ^^n-n=^^{nl + nl + nl) . A2 



Stating more precisely, in order for the space S to possess desirable properties, varieties 
of objects in the theory should be rich enough; viz. for a given object O there should always 
be another object O' which is an infinitesimal modification of O, so that the spectra for O and 
those for O' are infinitesimally close to each other. This condition is satisfied in the case of Sn', 
for any given geometry, there always exists its infinitesimally modified geometry Q. 
'See the beginning of §|| for more details on the notations such as #((?). 
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We recall that we have introduced the eigenfunction f^'^ and the eigenvalue 



Xa of the Laplacian A on the regular 3-torus T^(l) (see after Eq.(|l2|)). They are 
related to f^^ and as 

fP ■■= /f , := a'Xn . (A3) 

For notational simplicity, we often write /a and Ja instead of and f^, 
respectively, wherever explicit indication of the label (n, a) can be omitted. 

We recall the conformal equivalence between and T^(l) (see the argument 
after Eq.(|l^)). Because of this equivalence, we can define the integral over T^(l) 
for any function A and any symmetric tensor field Aab on S^, 

{A)~^ ■■= f UA U d'x , {Aab) A ■=!- f d'^U Aab d'U d^x . (A4) 
Then, it is useful to introduce the coefficients c'aa'^ '^aa' ^^'^ ^AA' '^^ 

CaA' ■ = {Ma = ^{fA')A 

- r_^#(<?ii)^r./2^#("')A .A.,,. A . 

[i^dadbfA')A =V^{^dadbfA')A (^'^0) , 

{n' = 0) 



'^AA' 



-2 (1) 



AA' 

C^AA' '■ = ~{fA'hab)A = -^{fA'hab) a 

= {l-2nl)c% . (A5) 

Here ny {n±) is, roughly speaking, the component of n parallel (perpendicular) 

to n' ; n| := (n||/|^|^ and := {n±/\n\Y, so that n'j and n]_ possess only 
the information on the relative direction between ft and n'. Thus it is obvious 
that the coefficients Caa'^ c^aa' ^'^d c^aa' do not depend on the information of |n|, 
\n'\ but on the information of the relative direction between n and n' . 

In particular, we note Caq = —Cao = 1 ^^id Caq = by definition. It is also 
helpful to note that, when n' 7^ 0, 



I ~. 7. 7 r^. I 7 

f .,h . \ ^ 



= {^^fA'hab) A = Va^i^AfA'hab) 
AA' ^A' 



There is an identity 



SI) ^(2) . (3) 



^AA' 



^c'a'a' + c'a'a' = , (A6) 



(3) 

so that one among the three kinds of coefficients, say Caa'i can always be elimi- 
nated from the formulas. 



See the beginning of §^ for notations. 
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A.l Scalar modes 

From the basic modes Ja in Eg . ([Al]) , the components of the perturbations jab in 
Eg . (11^) can be produced. 
Let us define dab as 



_ 1 „ „ ^ 1 T ^ r 

'ab ' 



^ab '■= i-dadbfA - ^AfAhab ■ 
AA 'J 



Then, fAhab and o"^ form the basis of the scalar modes of the metric perturbation 

lab- 

7r = I E 7^^^ fA Kb + T.'l''2 , (A7) 

A A 

where 7'^' and 7'^' are arbitrary expansion coefficients. 
We note some relations including a^, 

D^'ai = -\djA' , D''&~ai = '^-\A'fA' , 

^^ab = -^A'Ia' , 
K'^ablA — ^AA' r^^AA' — ^\^AA' ^ ^AA') i I 

where Eg . ( |A6|) has been used in the last line of the fourth formula. 

Then, it is straightforward to compute from fi^^^^'^^^^^ ■= a^{^3ab + 7!^'*''^'^'') the 
curvature guantities up to 0(7), 



^(scala.) ^ V'(a^A^)(7'^^ + 7'2) h hab - + I'Wbh 

^ A ^ A 

^(scalar) ^ + 1^) f A , 



^'a 

-r^^^ = -lY'i<^'^AW2+i''2)&^b ■ (A9) 

A. 2 Vector modes 

We choose the basic solutions of the eigenvalue problem, Ava = —Xva, with 
D'^Va = 0. We get the solution 

"^'^-^ = a fin,,A.ei.^M'\ , (AlO) 

where (e''-'^)^ := Sja can be identified with the unit vector in R'^ along x-'-axis 
(j = 1, 2, 3), and {a{e^^^)°'}j=i^2,3 forms the orthonormal bases in the t = constant 
geometry.0 They are normalized as 



Note that the indiees are lowered by hab in Eq.(|lJ) and raised by its inverse, h 
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Prom v^''^'^\ we can construct the traceless tensor C, 



ah 



[DaVb + 



DbVa), which is normahzed as {Cab'''^\ Ccd ^) = ^j,j'^ii,n'S<f,a' ■ 

Then, Cat'^'' :— Va-C^f'^^ form the basis of the vector modes of the metric 
perturbation jab- 

^(vecto.) ^ ^'y.) ^A, (All) 

A,j 

Wherever the detailed index structure is not essential, let us write A instead of 
Aj for notational neatness, like 7^', and I]'a7^a Cab- show that 



The curvature quantities up to 0(7) turn out to be 



(A12) 



(vector) 
■ab 



(A13) 



A. 3 Tensor modes 

We choose the basic solutions of the eigenvalue problem, Awab — —^Wabi with 
D^'Wab — 0, — 0. We get the solution 



w 



^ (ft. a, a) 



ab 



(a) 

ab ' 



(A14) 



where 



Kb — 



^ab 



1 
1 

71 



/ 



V 1 



/ -1 



(11) 



V2 
I 



V2 



, ^ab 



{III) 



\ 1 




and e<i) = = (1, 0, 0), e<2) = ^^^) = (0, 1, 0), e<3) = ^^^^) = (0, 0, 1). 
They are normalized as 



ab^cd^{n',a',a') 



bd 



^a,a' ^n,n' ^ff,S' 



Then, - ^-(^'"^ 



perturbation ^ab'- 



-^w"^' "' form the basis of the vector modes of the metric 



-(tensor) _ ST^ IT) ~Aa 
— 2^1 Aa ^ab 



Tab 



A,a 



Let us introduce the coefficient c^Tl,„ as 



AA'a 



A ■ 



(A16) 



(A17) 
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For a = i = 1,2,3, 



JT) 

(n,(7){m,(7')i 



(1 - '^n,,o)^^(^f+l - '^'i+2)(^{n,<j){m.,S') 

(when at most one of ni,n2,n^ vanishes) 


(when at least two of ni,n2,n^ vanish) 



(T) 
AA' 

The curvature quantities up to 0(7) become 



For a = I, II, III, it turns out to be c)^^/„ = 0. 



(A18) 



r?(tensor) _ ^ V^' x ^.(T) ~A D(tensor) _ (tensor) _ ^S'^KtA (^^Q\ 

A ^ A 

where we have used a simphfied index notation A for Aa for notational neatness. 
We follow this convention wherever the detailed index structure can be omitted, 
and we write simply, e.g., 7^', w^^ and Y.aIa^ "^ab- 



B Basic results of the linear structure-formation 
theory 

We summarize the time evolution of the linear perturbations for the typical cases 
when (I) Aphys ^ c^t and 6^ 7^ 0, (!') Aphys ^ c^t and 6* = and (II) Aphys <^ Cgt, 
where Aphys is the physical wave length of the perturbation mode [Q. 



B.l Behavior of perturbations: The case of Aphys ^ Cgt and 

A perturbation mode satisfying Aphys ^ Cgt does not oscillate sufficiently during 
the cosmological time-scale t. It is a long-wave length, slowly varying fluctuation 
compared to t. Let r be the oscillation time-scale of the perturbation, r = 
Aphys/cs. Then t -C r, so that the only important time-scale in this situation is 
t, and perturbations naturally vary as a power of t. 

We assume v = constant = for simplicity. 

In this case, 

, , 2(1 + 3z/) , 
S{x,t) = F(f) ^ 7(f,t) = J- -{S , 



u{9u - 1) 

Z/)(6Z/ +1)5 1 2(2-3'^) 
OC — t 

3u(l + uY t V 



9{x,t) = ^^TTTT-^l oc-rww . (Bl) 



Here F{x) is an arbitrary function determined by the initial condition imposed 
at a certain time t = toi the factor l/u has been shown in the last term of the 
last line to indicate the singular behavior for u I 0. 
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As a typical situation, it is illustrative to consider the plane-wave type com- 
ponent of perturbation which is moving in the x^-direction. In this case, the 
behavior of 733 and 711 ~ 722 are 

8(6z/+l) a" ^ l/AphysV. 1-1^ 

where k is the coordinate wave-number {k = 27ra/Aphys)- 

It is convenient to summarize the relation of the order of magnitude of several 
quantities in a symbolical manner as 

(1^1 ^ \5\ ~ \t9\) ~ (|t7| ~ \t5\ ~ \t'e\) « I7./I ~ \t%i\ . (B3) 

Regarding the rate of change of 711, we see that t|7ii| ~ |7ii|- It means that 
\jki\ decreases prominently within the cosmological time-scale. It implies that 

l^i^i^'fl ~ {^^) I'Tiil^m)!" l7ii(^m)l< l7ii(^o)| where is the time when the 
perturbation of scale Aphys "enters within the horizon" due to the enlargement 
of the horizon scale (~ ct), viz. Aphys ~ Cgt. It indicates the linear perturbation 



theory is valid even after t = tin |14 



B.2 Behavior of perturbations: The case of Aphys ^ Cgt and 

^ = 

In this case, 

6{x,t) = F{x) tTOT , 7(x,t) = -Y~^^ ■ ^^^^ 

For the plane-wave type component of perturbation moving in the x^-direction, 
the behavior of 733 and 711 ~ 722 becomes 



_ 8(9// + 5) a" /AphysA . 1 ^0 



(B5) 



On the other hand, it turns out to be 



We can summarize the relation of the order of magnitude of several quantities 
symbolically as 

(I7I ~ \5\) ~ (inl ~ \tlui\ ~ m « \lki\ , (B7) 

The peculiar feature of the case 9 = compared to the case 6* 7^ is that \'^ki\ 
does not change prominently within the cosmological time-scale (t|7ii| ~ \5\ <^ 
l7ii|)- 
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B.3 Behavior of perturbations: The case of Aphys <C Cst 

A perturbation mode satisfying Aphys ^ Cst oscillates sufficiently during the cos- 
mological time-scale t. It is a short-wave length, fast varying fluctuation com- 
pared to t. Now t S> r (r = = ^ characterizes the oscillation period), so 
that we naturally pay attention to the process whose time-scale is 0{r). 
In this case, various quantities oscillate with the time-scale r: 

5 ^ -7 --733 0c(^^^J exp-.y -^dt , 

e ~ — —s , 

I -I II I I II f P + P \^^^ 1 1-3./ 1 l-Sf 

\5\ ~ 7 ~ 733 ~ ''■ F oc 7^ — (X ^=a 2 oc —;=t , (B8) 

where we have assumed p ~ Pav oc a"^^^"'"'^) oc t~^. 

We can summarize the relation of the order of magnitude of several quantities 
in a symbolical manner as ("~" indicates that the difference between the both- 
hand sides is at most o{6).) 

Aphys 
Cgt 



~ |7|~|733|)~(t|5|~t|^|)<^<1 
I711I ~ I722I ~ (o((5)) , 



^i7..i) -[^] \^\« \^\ ■ m 



2 



Cgt 
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